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Abstract 
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1 Introduction 



Polyakov [T] introduced monopoles in 2+1 dimensions as instanton tunneling events in com- 
pact gauge theories. The prohferation of these monopoles leads to confinement and to the 
absence of a Coulomb phase in such gauge theories, provided there are no gapless matter 
fields that can suppress the monopoles. In condensed matter physics, two-dimensional lat- 
tice quantum antiferromagnets can be written as compact U{1) gauge theories at strong 
coupling [2]: here, monopole events are accompanied by Berry phases [3}|4], which are re- 
sponsible for valence bond solid order in the confining phase jijjs] . 

We can also consider monopole operators at conformal fixed points of 2 + 1-dimensional 
gauge theories [6 -13 . These are gauge-invariant primary operators that determine impor- 



tant aspects of the structure of the conformal field theory (CFT). In the application to 
antiferromagnets, the scaling dimension of the monopole operator determines the power-law 
decay of the valence bond solid order at "deconfined" quantum critical points [9| [l0|[l4[[l5] . 

This paper will consider a different class of monopoles in 2+1 dimensions. We consider 
CFTs with a global U{1) symmetry. The CFT may also have fluctuating gauge flelds, but 
these play no role in the construction of such monopoles. Instead, the monopole is introduced 
by a background U{1) gauge fleld that couples to the CFT conserved current. A monopole 
with charge q inserted at r = vq, which we henceforth denote by A^g(ro), corresponds to a 
background gauge field configuration whose field strengthj^ 

ffiu = d^ay - dyU^ (1) 

integrates to 27rg over any small two-sphere surrounding the insertion point: 

/ = 27rg. (2) 

As we will see explicitly in Section [2| each such background monopole comes associated with 
a Dirac string that starts at r = tq. If the matter fields have integer U{1) charges, the Dirac 
string is not observable provided that q is an integer. 



Such monopoles appear to not have been considered until recently 16, 17 . They do not 
correspond to operators in the CFT in a strict sense; instead, they should be rather thought 
of as non-local background sources to which we couple our CFT. Studying the response of the 
CFT to such background sources provides useful information about the CFT, which can be 
used, for instance, to test various dualities |16j. In addition, these monopole insertions have 
been argued to play a crucial role in the structure of the compressible quantum phases that 

^In standard vector notation, instead of / we would use the magnetic field /3 — */, which can be also 
written as /3 = V x a. Eq. ^ becomes J^a • dA — 27rg, where dA is the oriented area element. 
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are obtained when a non-zero chemical potential is applied to the global U{1) charge 17, 18 . 
Specifically, they serve to quantize the U{1) charge and to determine the lattice spacing of 
Wigner crystal states such that there are an integer number of particles per unit cell; they are 



also important in determining the period of Friedel oscillations 18 , 19 of compressible states 
that do not break translational symmetries and may have "hidden" Fermi surfaces 17 20 -22 . 
In this paper we will restrict our attention to CFTs to which no external chemical potential 
has been applied. 

The definition oi Aiq presented above is imprecise, partly because the condition ^ 
does not specify / uniquely, and partly because we have not specified the allowed behavior 
of the charged matter fields close to the singularity at r = tq. Just as in the case of 
monopole operators in gauge theories [T], a precise definition can be given through the state- 
operator correspondence, or, more precisely, through an extension thereof to the present 
case. According to the state-operator correspondence, any local operator of a CFT inserted 
at the origin of corresponds to a normalizable state of the CFT on S"^ x M, where the 
M coordinate is interpreted as Euclidean time. A monopole insertion A4q is by no means a 
local operator, but it can nevertheless be defined as corresponding to the vacuum on S'^ (as 
opposed to any other excited state) in the presence of q units of background magnetic fiux 
(as in (|2])) that is uniformly distributed throughout the S'^. 

The monopole insertion defined above is a Lorentz scalar. It also has a well-defined scaling 
dimension Ag in the following sense. If we consider a background gauge field configuration 

corresponding to a monopole of strength g at r = ri and one of strength —q at r = r2, 
the partition function in the presence of these two monopole insertions has power-law decay 
with the relative distance |ri — r2|, namely 



, „ fP^exp (- r d^xCla]) 



°^ ] ' (3) 



where we denoted by C[a] the Lagrangian of the CFT coupled to a^. One can extract A, 
from the exponent in ([s]). Equivalently, in view of the definition of Aig through the state- 
operator correspondence described above, one can also map a single monopole insertion on 
M.^ to S*^ X M and identify A^ with the ground state energy on 5*^ in the presence of q units 
of background magnetic fiux. Explicitly, we have 

A, = = -logZ,, (4) 

where Zq is the partition function on S"^ x M, and J-'q the corresponding free energy. 

The question that we will address in this paper concerns the scaling dimensions Aq of 
the monopole insertions Aig in simple CFTs. To calculate these scaling dimensions, we will 
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use Q. For certain free CFTs with global U{1) symmetry, one can infer Ag from existing 
results in the literature. A simple example is the free CFT of Nj complex fermions. The 
Lagrangian 

Nf 

Cf = J2^i{^0)^l^a (5) 

a=l 

is invariant under a U{Nf) global symmetry under which ipa transforms as a fundamental 
vector with Nj components. We can consider the diagonal U{1) subgroup of U{Nf), which 
we couple to a U{1) background gauge field such that the modified Lagrangian is 

Nf 

Cf[a] = J2i^i{z$ + <^)ij,. (6) 

a=l 

As above, we consider monopole insertions of q units of background magnetic flux. Using 
Q, Ag can be computed from the partition function on S'^ x M, which is now a Gaussian 
integral because the Lagrangian ^ is quadratic in %jja and there are no interactions. The 
same Gaussian integral was calculated in Ref. [T] as part of a slightly different problem: 
The authors of Ref. [7] were interested in computing the scaling dimensions of monopole 
operators in three-dimensional QED with Nf flavors, which is the same theory as ([6]), with 
the exception that the gauge field a^j would be dynamical. While the leading large Nf result 
of Ref. f?) is only approximate for QED (because there are corrections coming from the 
fluctuations of the gauge field), in the free fermion theory ^ one obtains an exact result 
that holds at all Nf^ We reproduce the dimensions for the first few lowest values of q 
in Table [T| Similar results for a non-supersymmetric free theory of Ni, complex scalars are 



Q 


^JNf 








1 


0.265 


2 


0.673 


3 


1.186 


4 


1.786 


5 


2.462 



Table 1: The scaling dimensions of the monopole insertions A4q in the free theory of Nf 
fermions corresponding to the diagonal U{1) subgroup of the global U{Nf) symmetry group. 
These results are exact in Nf. 

presented in Appendix |Aj 

^We should restrict to Nf even in order to avoid a parity anomaly. 



3 



In this paper we are interested in the more comphcated case of an interacting CFT 
with global U{1) symmetry. The simplest such CFT is the XY model, described by the 
Wilson-Fisher fixed point of the 0"^ field theory of a complex scalar field (p. Starting with 
the Lagrangian 



-XY 



\d^^\^ + s\<j>\' + u\<j>\\ (7) 



the Wilson-Fisher fixed point is reached in the infrared provided that the coefficient s is 
tuned to zero. This theory has a global U (1) symmetry under which is rotated by a phase, 
and it is this U{1) symmetry that we couple to a background gauge field a^. The Lagrangian 
in the presence of is 

CxY[a] = 1(5^ - ta^W + + ■ (§) 

As in the previous examples, we can consider a monopole configuration with q units of 
background magnetic fiux as defining the insertion Aiq. 

Unfortunately, the Wilson-Fisher fixed point of a single complex scalar cannot be accessed 
perturbatively, so we will compute the dimensions by first generalizing Cxy [«] to a theory 
with Nf, complex scalars with Lagrangian 

C[a] = Y,\i^,-^aM' + s\M' + u{\^\') , |0|2 ^ ^ |0,|\ (9) 

a=l a=l 

and then performing a 1/Nh expansion. Our goal in this paper is to find the first two terms 
in this expansion. The CFT (obtained by setting a = in ([9]) and tuning s to zero) is the 
Wilson-Fisher fixed point with 0{2Nb) symmetry. The U{1) symmetry that we consider is 
a subgroup of 0{2Nb) that acts by rotating each complex scalar by the same phase. 

The rest of this paper is organized as follows. In Section [2] we set up our conventions and 
explain the method we use to compute Ag in the model ^ in more detail. In Section |3] we 
perform the leading order calculation in N^. To this order, we find agreement with the results 
of Refs. [6, 11 on the leading large Nb dependence of the dimensions of monopole operators 



in the CP^'~^ model. Indeed, to leading order in Nb, one can ignore the contribution to the 
S"^ ground state energy coming from the gauge field fiuctuations in the CP^*"^ model, so 
the scale dimensions of the monopole operators in that model should agree with those in the 
ungauged theory (joj). In Section [4] we compute the leading 1/Nb corrections to Ag. We end 
with a discussion of our results in Section [H 
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2 Method 



We consider the 0{2Nb) scalar field theory defined on an arbitrary conformally fiat manifold 
by the action 



a=l 



n 



(10) 



where TZ is the Ricci scalar of the background metric g^y, and is a background gauge 
field. The only dynamical fields are the complex scalars 0^ and the Lagrange multiplier field 
A. It can be checked explicitly that this action is invariant under the Weyl transformations 



0a ^ /(r) 



-1/2, 



111 



for which / can be taken to be an arbitrary real-valued function. We will be interested in 
the action (10) on two conformally fiat backgrounds: and 5*^ x M, which have 7^ = and 
TZ = 2, respectively. 

On ]R^ in the case where and = 0, the action (|10| describes the Wilson-Fisher fixed 



point of 2Nh real scalars. Indeed, one can add the term J d^rX'^/{4u) to this action without 
changing the IR fixed point, because u fiows to infinity; integrating out A produces the 
interacting theory (|9| with = s = 0, which represents the more conventional description 
of the Wilson-Fisher fixed point. The monopole background (a^) = AJ^ that corresponds to 
an insertion of Aig at the origin of satisfies, in spherical coordinates^ 



dA" = I sin A 



;i2) 



which follows from We can work in a gauge where 



A'^ir) = |(1 - cos 



(13) 



This background gauge field is well-defined everywhere away from 9 = tt where there is a 
Dirac string. This Dirac string is not observable provided that q is taken to be an integer. 
Starting with the theory on M? in the monopole background (13), the theory on S"^ x M 



can be obtained from a Weyl transformation as in (11). Indeed, writing the fiat metric on 



In standard vector notation, we would write V x yt"? = qerl(2 \r\ ) instead of Q, and A'' = |(1 - 



cos0)/(rsin6')e0 instead of (13) in flat space. On S' 



I, we have ^ ^{l - cos 6')/(sin 6')e^ 
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in spherical coordinates as 



ds^ = dr^ + r^dn^ , (14) 
and defining r = e'^, we obtain a metric conformal to S"^ x M: 



ds^ = {dr^ + dn^) . (15) 

So if we send g^l — > e'^'^g^l = g^l^^ and at tlie same time rescale (pa e'^/^^a, A — )■ e^^A, 
— > as dictated by (11), we obtain tlie action on S"^ x M. Tlie monopole background 



(13) now corresponds to a constant magnetic field uniformly distributed over S"^. 

As explained in the introduction, we identify the scaling dimensions Ag of the monopole 
insertions A^, with the ground state energy J'g on S"^. We expand this ground state energy 
at large A'';, as follows: 



= NtT^ + + 0{1/Nt) . (16) 

When g = the operator Aig is just the identity operator and it corresponds to the ground 
state on S"^ in the absence of any magnetic flux. We expect this operator to have vanishing 
scaling dimension. Indeed, we will check explicitly that J-'q = in our regularization scheme. 

We now turn to the evaluation of J-'^ in the next section and of SJ-'q in Section |4j We 
will work solely on S"^ x M whose coordinates we denote collectively by r = (r, 6, 0). 



3 Ni) = oo theory 

In computing the leading large A^";, contribution to the ground state energy on 5*^, one can 



evaluate the partition function corresponding to (10) in the saddle point approximation 
where the fluctuations of the Lagrange multiplier field A can be ignored. However, A should 
be adjusted such that the ground state energy is minimized. We thus expand the Lagrange 
multiplier about its saddle point value a^ 

tX = al + ^+t~X, (17) 

where will be determined shortly by the saddle-point condition, and A is a fluctuation 
that we will consider in the next section. 

^This notation has been chosen to be compatible with Ref. jllj that studied the CP^*"^^ model. 
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We expand the field 0^ in terms of tlie monopole liarmonics defined in Ref. 23 :^ 



l=q/2 m 



27r 



Zem 



q/2,im\'7, • 



;i8) 



The quadratic action for the (pa then takes the diagonal form 



Ni oo £ 
a=l e=q/2 m=-e 



2tt 



[u^ + {e+l/2f + al] \ZemA^)\' 



(19) 



where we have used the fact that the eigenvalues of the gauge- covariant Laplacian on S'^ are 
£{i + 1) — {q/2Y |23]. From (19), it is easy to read off the leading approximation to the 
ground state energy at large A^";,, which comes from performing the Gaussian integral over 
the scalar fields 0a, or equivalently over the coefficients Zf„i,a- The coefficient J-"^ appearing 



in (16) is then 11 



= / ^ E (2^ + 1) l°g y + + 1/2)^ + . 



(20) 



This expression is divergent, but it can be evaluated, for instance, using zeta function regu- 
larization. First we write formally log A = —dA~^/ds\^_^ in all the terms of (20), then we 
evaluate the sum and integral at values of s where they are absolutely convergent, and at 
the end we set s = 0. Performing the u integral, we obtain 



(21) 



s=0 



which still diverges when evaluated at s = 0. We then use the identity 



i=q/2 



(£ + 1/2)2(1-^)+ (^-s)a^(£ + 1/2)2^ 



s=0 



g(i - (f 

12 



2 



(22) 



This identity can be derived by writing the sums on the left-hand side in terms of the Hurwitz 
zeta function CiSiO) = Yl'^=o^/i''^ + '^)* ^^"^ analytically continuing to s = 0. The terms 
on the left-hand side of (22) are nothing but the large i expansion of the terms in (21), so 
subtracting (22) from (21) yields a finite result when s = 0. Adding and subtracting (22) 

'^Note that our definition of q difl^ers from that of Ref. 23 by a factor of two. 
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from (21), we tlierefore find 



oo 



i=q/2 >- 



(£ + 1/2) [(£ + 1/2)2 ^ ^2] 1/2 _ ^ ^/2) 

g(g2-l) ^ 



2 - -a^ 



(23) 



24 



+ 



wliicli involves a convergent sum over C. that can easily be evaluated numerically. 

The value of is not arbitrary, but should be chosen so that the saddle point condition 



(24) 



is satisfied. In our case, where J-'?° is given by (23), we therefore have 



E 



+ 1/2 



2 



(25) 



1^2 VV'(^+ 1/2)2 + a2 
For the first few small values of q, we give in Table [2] the solutions of this equation as well 



as the corresponding values of J-'^ obtained after plugging these solutions back into (23). 
The values of J-'^ agree precisely with those obtained in Ref. joj in the large Ni, limit of the 



Q 




■TOO 











1 


-0.4498063 


0.1245922 


2 


-1.3978298 


0.3110952 


3 


-2.8454565 


0.5440693 


4 


-4.7929356 


0.8157878 


5 


-7.2403441 


1.1214167 



Table 2: The values of that solve (25) and the corresponding coefficients J-'?° that enter 



the large Nb expansion of the ground state energy (16) on S . 



case of these results is notable. The value Og 



The q 

the conformal mapping between 
at intermediate stages. 



and 5*2 x 



Also Jo°° 



= is just that expected from 
0, a result that was not evident 
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4 l/A/^ corrections 



4.1 General structure 

The leading 1/N}, correction to the result of the previous sections comes from the contribution 
to the S"^ ground state energy coming from the fluctuations A of the Lagrange multiplier. 
Let us begin by discussing the general structure of this correction. 

After integrating out 0^, the effective action for the fluctuations takes the form 



\ I rSr' ^^^^\{r)D%r, r')A(r') + . . . , (26) 



where we omitted higher order terms in A. The kernel D'^{r,r') appearing in eq. (26) is 
nothing but the two-point correlator of \4>af '■ 

D\r, r') = (|</>„(r) I' |</),(r') |') = iV,G(r, r')G*(r, r') , (27) 

where we introduced the Green's function G{r,r') = {(f)*{r)(f){r')) for a single complex field 
(j) in the background monopole flux We will compute this Green's function shortly. 



Because of the explicit factor of A*";, in (27), at large A*";, we can ignore the higher order terms 
in (26), and evaluate the contribution from A to the partition function in the saddle point 
approximation. The coefficient 6J^q appearing in eq. ( 16 ) can then be obtained by performing 
a Gaussian integral, which yields 



5jr = llogdetD". (2^ 



To calculate logdetD'^, we should diagonalize the kernel D'^. This diagonalization is 
accomplished by expanding A and D"^ in terms of the appropriate spherical harmonics. These 
quantities do not experience a net monopole flux, because they are neutral, and so we 
(fortunately) do not need the monopole spherical harmonics here. The expansions 

A(r)= f ^e--r,„(e,0)A,„(u;), 

r Z (29) 

Im 

yield a diagonal effective action 



Scs — - 



\ I ^E^'H|A.™(u;)p. (30) 
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Eq. (28 



) then gives 



5Fc 



1 f du 



(31) 



In the following subsections we present expressions for the kernel in (30): We first present 



the simpler kernel at g = 0, and then the kernels at general q. We will check explicitly that 
= 0, as required by conformal invariance in the absence of any monopole insertions. 



4.2 The kernel of fluctuations at g = 

At g = 0, it is not hard to obtain the Green's function on 5*^ x M starting from the Green's 
function on M^, l/(47r|r — r*|), and using the conformal mapping explained around equation 



(15). The result is 



G(r,r') 



47r Y/2(cosh(r — r') — cos 7) 
where 7 is the relative angle between the two points on S"^ defined through 



(32) 



cos 7 = cos Q cos Q' + sin Q sin Q' cos(0 — 0') 



(33) 



Using (27), (29), and (32), we obtain 



N 1 r , r ■ n,n e^'^"P/(cOS( 



IGtt 



00 Jo 



(cosh r — cos I 



(34) 



We performed these integrals analytically for a number of small values of £; from the structure 
of these answers we deduced the general result: 



D^(tu) 



T{{i + l + tco) /2) 



4T{{i + 2 + iio) /2) 



which can be written more explicitly as 



tanh(7rC(;/2) 



uj coth(7ra;/2) 
8(a;2 + 1) 



fr (a;^ + (2n-l)' 

n=l ^ ' 

(w^ + An^) 



n 



L (u;2 + (2n + 1 



\2\ ■ 



(35) 



(36) 



In the limit of large u and i we expect the A self-energy to be given by the flat space 
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limit 



1 



1 



Stt^ p"^{p + p')'^ 8p 



Indeed, expanding (35) with the help of the Stirling approximation, we find 



i{i + 1) 



+ 



8y/uj^ + £{£+!) 32 + + l)f^ ' ^ \^ {co^ + £{£ + ^ 
which agrees with (37) upon using the identification p ~ ^u^ + £{£+!). 



(37) 



(38) 



4.3 The kernel of fluctuations for general q 

Now we turn to the much harder case of non-vanishing q. In this case we don't have a simple 
closed form expression for the scalar Green's function, so we turn to the mode expansion 



(18). From the action (19) we deduce that the Green's function for a single (pa is 



oo „ 

G(r,r')= J2 / 



— f 
27r 



,ja;(r— r') 



En 



q/2/rn{^^4>)yq/2,lra{0'A') 



_m=-i 
-Ege\r-T'\ 



u;2 + (£ + 1/2)2 + a2 



(39) 



2E, 



ql 



In writing the second line we defined the energy 



',, = ^(£+l/2)2 + a 



2 

9 ' 



(40) 



and performed the uj integral; we also performed the sum over m, which, up to the phase 
factor 



J0 



1 



cos(^/2) cos(^72) + e-^('^-'^') sin(^/2) sin(^72) 



(41) 



cos(7/2) 

discussed in Ref. [24], yields a polynomial in cos 7 that can also be written in terms of the 
monopole harmonics as 



FqAl) 



2£+l 
An 



Yq/2A-q/2i7,0) . 



(42) 



(See Appendix |B] for more explicit expressions for Fq^£('y).) Here, 7 is the relative angle of 
the two points on 5*^ defined in (33). 



From (27) and (39), we can now determine D'^{r,r'). Further extracting D'^{u)) using 
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(29) we obtain 



Dj{uj){27T)6{uj + uj') 



{21 



1 /. 

'pi III— n 10 



l.'/"=q/2 ■ 

^-{Eqtl +E^ill )\T-T'\-iWT-iuj'T' 

X 



We can simplify this expression to 



{2e + 



I) ^ 

' eei=q/2 



^EqlEqli 



Eqii + Eq£ii 



(43) 



2Eq£iEqeii{u'^ + {Eqii + EqiiiY) 



lDii,i',i" 



where 



Jo 



(44) 



(45) 



This is an integral of three monopole harmonics and can be expressed in terms of the Wigner 
3-j symbols as 



(2£ + l)(2f + l)(2r + 1) 
32^^3 



ffi nil 



-g/2 g/2 



(46) 



We can check, for instance, that this result equals (34) for g = and £ = 

tanh(7ra;/2) 



1 1 



(47) 



4.4 Numerics 



The results of the previous sections are all we need for calculating numerically the correction 



to the scaling dimensions of the monopole operators. Unfortunately, the expression (31 ) 



is formally divergent, as can be seen for instance in the case g = where we know D^ibS) 



explicitly, and hence eq. (31) is not suitable for numerical evaluation in its current form. 



However, we expect the divergences to be independent of g, so the differences ^Tq^ — ^J^q^ 
should be finite and shouldn't require regularization. Moreover, it must be true that bT^ = 0, 
because the case q = corresponds to an insertion of the identity operator, which should 
have vanishing scaling dimension. (See Appendix [C| for an explicit check that 5J-'o = 0.) 



Subtracting 5J-'o from (31), we can then also write 5J-'q as 



(45 



12 



which we evaluate numerically for the first few lowest values of q. 



In evaluating (48), one has to perform three sums (two when calculating D\{u:) using 



(44) and one in (48)) and one integral over u. Let us first comment on the two sums in 



(44). For fixed the sum over C" in (44) has only finitely many non-zero terms because the 
3-j symbols in (46) vanish unless £, I' and £" satisfy the triangle inequality. To see whether 
or not the remaining sum over I' is convergent, one should find an asymptotic expansion 
at large for the terms in this sum. While for general £ it may seem hard to do so, it is 
easier to first fix £ to a small value for which the sum over ^" has 2£ + 1 terms that can be 
written down explicitly, and the large ^' asymptotics can be easily computed. Repeating this 
procedure for several values of one can infer the large ^' asymptotics for all I by noticing 
that all the expressions involved are polynomials in ^{[ + 1). The first few terms are 



1 1 3 - + 2^{^ +l)-uj 



- ' -V- ' V (49) 

This expression shows that the sum over i' is absolutely convergent. To save computational 
resources, one can use a mix of numerical and analytical techniques in evaluating Dl{uj): the 
terms with low should be summed up explicitly, while for the terms with large I' one can 



sum up analytically the approximate expression ( 49 ) developed to a higher order of accuracy. 



(In our computations, we developed the large approximation up to order 



Lastly, in calculating (48) one should be wary that there could still be divergences. We 



find that imposing a relativistic cutofl^ 

a;2 + £(£+l) <L(L + 1), (50) 
yields a finite answer as we take L — )■ oo. The absence of divergences relies heavily not only 



on the choice of cutoff (50), but also on choosing the value of that solves eq. (25); for 
other values of there would be divergences. See Figure [l] for a plot of 5J^q in terms of 1/L 
in the case q = 1, where from the large L extrapolation we obtain 5J^i ~ —0.057. In this 
case we therefore conclude that the scaling dimension of the monopole operator A^i is 

Ai = 0.125A/'6- 0.057 + 0(l/iV5), (51) 

where we included the leading large Nf, behavior that was also given in Table |2} Repeating 
this procedure for the first few small values of g, we obtain the results in Table |3| This is 
the main result of this paper. 



^At high energies the Lorentzian theory has 50(2, 1) symmetry that is also obeyed by the cutoff (50), so 
the speed of Hght is not renormahzed. If one chooses a cutoff that breaks the 50(2, 1) symmetry, then there 
are finite corrections to the speed of hght in the IR that have to be accounted for. 
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5F 



-0.0571 



-0.0572 



-0.0573 - 



-0.0574 - 



-0.0575 




0.02 



0.04 



0.06 



0.08 



0.10 



1/L 



Figure 1: The coefficient 6J-'i evaluated numerically from (48) using the relativistic cutoff 



(50) as a function of the inverse cutoff scale 1/L. The solid line is a quadratic fit from which 
we extract the value SJ^i ~ —0.057 as we take L — oo. 





A, = 










1 


0.125iVfe 


- 0.057 + 0{1/Nb) 


2 


O.SllNb 


~ 0.152 + 0{1/Nb) 


3 


0.5UNb 


- 0.272 + Oll/Nb) 


4 


O.SieATfe 


-O.AU + oll/Nb) 


5 


l.UlNb 


- 0.575 + Oll/Nb) 



Table 3: The scaling dimensions of the first few monopole operators A^g in the Wilson-Fisher 



CFT of A*";, complex scalars in the large Nb expansion (16). The leading large A*";, behavior 
was computed in Section |3| and agrees with resuhs from the CP^""^ model |6j. The 0(A/'°) 
term was computed numerically using ( 48 ) . 



5 Discussion 



Following recent work 16, 17 , in this paper we considered monopole insertions in 2 + 1- 
dimensional CFTs that have a global U{1) symmetry. A simple example of such a CFT 
is the Wilson-Fisher fixed point of the XY model. Critical exponents of this CFT have 
long been the focus of much study, and are among the most accurately known non-trivial 
exponents of higher dimensional CFTs [25]. Associated with the monopole insertions, we 
have a new set of critical exponents of this venerable CFT. We computed these exponents 
[i.e. monopole scaling dimensions) to next-to-leading order in the 1/Nb expansion of a theory 
with Nb complex bosons. Our results for the scaling dimensions are summarized in Table [3] 
The numerical series in Table |3] appear to be reasonable even when evaluated at A";, = 1. 
It would be interesting to also compute the monopole scaling dimensions in Monte Carlo 
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simulations or series expansions, such as those in Ref. 
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A Free scalar theory 

We can also calculate the scaling dimensions Ag in the free theory of A*";, complex scalars. 
The only difference from the Wilson-Fisher CFT is that the action for the free theory does 
not have a Lagrange multiplier A, but there is a conformal coupling TZ |0|^ in the action, as 



in (10). The ground state energy on in the presence of q units of magnetic flux that we 



obtain by integrating out the scalars is N^J^^, where J-"^ can be computed from (|23|) with 



ttg = — as appropriate for conformally coupled scalars. See Table |4] for a few particular 
cases. These results are exact. 













1 


0.097 


2 


0.226 


3 


0.384 


4 


0.567 


5 


0.770 



Table 4: The first few scaling dimensions Ag of the monopole insertions A4q in the free CFT 
of Nh scalars. 



B Monopole harmonics 



We start with the relation 



m=-e 



Yg*/,^,je,<P)Yg/,Me',<p') = FgMe 



(52) 
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where F is defined in eq. (42) and the angles 7 and 6 are defined in eq. (33). 



Above, we have used the functions 



(2£ + i; 



(9,0) 



2-/^(^) (l + cos^r/^P,°:^^/,(cos, 



2-5/2 f?^](l + cos ^^"^/^-i 



(i + q/2)PlXn{^ose) + (£ - q/2 + I) PI^-'{cos 9) 



(£ + 1/2) 



(53) 



The special values are 
'2£ + l 

FqAO) 



,Pi{cose) ifg = 0, 

47r J 

71 ^°^^)"'^' [P,-1/2(C0S^^) + P,+i/2(cos^^)] if g = 1 , 



(54) 



etc. 



C Calculation of 5J^o 

We now show that using zeta-function regularization we find 6J-'o = 0. Using the infinite 



product representation for the hyperbolic tangent and cotangent in (36), one can show that 

00 

log D°^{uj) = ^ (-l)''+^log(w^ + k'^) + (^-independent terms) . (55) 



k=e+i 



The cj-independent terms do not contribute to 6J^o in our regularization scheme. With the 
help of 



du 
2^ 



log(w2 + a^^ 



(56) 



which can be derived, for instance, by rewriting (|56|) as 
d r du 1 



ds J 2tT (w2 + q2) 



d y/7TT{s - 1/2) |a 



l-2s 



ds 



s=0 



27rr(s) 



(57) 



s=0 
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we can perform the u integral in (31) and we obtain 



5:F, = -Y.^~lY{2^ + ^) Y.{-ifk. (ss) 

1=0 k=l+l 

The sum over k can be written in terms of the Hurwitz zeta function (^(s, a) = Xl^o / {f^+c^Y 
as 

oo 
k=l+l 

SO then 
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